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The (3+l)-dimensional ideal hydrodynamics is used to simulate collisions of gold 
nuclei with bombarding energies from 1 to 160 GeV per nucleon. The initial state 



O 

Q 

' is represented by two cold Lorentz-boosted nuclei. Two equations of state: with 



and without the deconfinement phase transition are used. We have investigated dy- 
namical trajectories of compressed baryon-rich matter as functions of various ther- 



^1^' modynamical variables. The parameters of collective flow and hadronic spectra are 

. calculated. It is shown that presence of the deconfinement phase transition leads to 



increase of the elliptic flow and to flattening of proton rapidity distributions. 



I. INTRODUCTION 

Relativistic hydrodynamics is a very popular approach to describe high-energy nuclear 
colhsions. Especially useful is its sensitivity to the equation of state (EoS) of strongly 
interacting matter and, in particular, to its phase transitions. In fact, extracting this EoS 
is the main goal of heavy-ion experiments. Some signatures of a deconfined quark-gluon 
plasma (QGP) have been already found in RHIC experiments with cm. bombarding energies 
yJSNN = 60 — 200 GeV. Presumably, this plasma can be also created at lower (SPS, AGS) 
energies. A more detailed data should be obtained in the low-energy runs at RHIC and 
future FAIR and NICA experiments. 

A large amount of experimental data on nuclear collisions at AGS, SPS and RHIC energies 
have been successfully described by different versions of the hydrodynamic model. The 



first model of this kind has been proposed by Landau One can roughly divide the 
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existing versions of ideal hydrodynamics into two classes. The first class includes the models 
which apply fluid dynamical simulations from the very beginning, i.e. starting from cold 
equilibrium nuclei. The models of the second class introduce an excited and compressed 
initial state - a locally equilibrated "fireball" formed at an early non-equilibrium stage of the 
collision. Up to now many versions of the fireball-based hydrodynamic model were developed 
ranging from simplified (1+1)- and (2+l)-dimensional models to more sophisticated (3+1)- 
dimensional ones (see u\ and references therein). 

Historically, early 3D hydrodynamic models of relativistic nuclear collisions p|, |^ used 
cold Lorentz-contracted nuclei in the initial state. However, such models become less justified 
at high bombarding energy due to the nuclear transparency effects. 

One should have in mind that the hydrodynamical approach can not be directly applied 
to late stages of a heavy-ion reaction when binary collisions of particles become too rare to 
maintain the local thermodynamic equilibrium. A usual way [5J to circumvent this difficulty 
to introduce a so called "freeze-out" criterium to stop the hydrodynamic description at 
late times. A more consistent procedure has been proposed in Refs. 0, Q] within a hybrid 
"hydro-cascade" model. 

In this paper we formulate a version of the ideal (3+l)-dimensional hydrodynamics suit- 
able for the domain of NICA-FAIR energies. This model belongs to the first class and uses 
a new EoS 8| with the deconfinement and liquid-gas phase transitions. By performing sim- 
ulations with and without the deconfinement phase transition we try to find the observables 
which are sensitive to creation of the QGP. A more detailed description of the results will 

n 

be given in Ref. |9[. 

II. FORMULATION OF THE MODEL 

Below we study the evolution of highly excited, and possibly deconfined, strongly- 
interacting matter produced in ultrarelativistic heavy-ion collisions. It is assumed that 
this evolution can be described by the equations of ideal relativistic hydrodynamics jlo| : 

d,T^'' = 0, (1) 
d^{nu^) = 0. (2) 
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Here T^^" is energy- momentum tensor {h = c = 1) 

T^"" = {e + P)u''u'' - Pg^"", (3) 

e, n and P are the rest-frame energy density, the net baryon density and pressure of the 
fluid, = (7, 7t>)'^ is its collective 4-velocity, 7 = (1 — v is the fluid 3-velocity, g^'' 
is the metric tensor. To solve Eqs. (l)-(3) it is necessary to specify the EoS P = P{n,6) of 
the fluid and initial conditions. 

In our calculations we use the EoS of strongly interacting matter with the first order 
deconfinement phase transition (EoS-PT) . The hadronic phase is described as the hadron 
resonance gas including known hadrons with masses up to 2 GeV. Finite size effects are 
taken into account by introducing the excluded volume corrections. The same excluded 
volume parameter Ve = 1 fm^ is used for all hadronic species. A Skyrme-like mean field 
U = U{n) is added to guarantee that the hadronic matter has correct saturation point and 
the hquid-gas phase transition. The quark-gluon phase is described by the bag model with 
lowest-order perturbative corrections. The phase transition boundaries and characteristics 
of the mixed phase (MP) are found by using the Gibbs conditions with the strangeness 
neutrality constraint. 

To probe sensitivity to the EoS, we have also performed calculations with the EoS of 
ideal hadron gas (EoS-HG). In this case we disregard the excluded volume effects assuming 
that fe = 0. To provide stability of initial nuclei we also introduce the Skyrme-like mean 
field. The parameters of EoS-PT and EoS-HG are given in Ref. 8|]. Our analysis shows 
that although generally the EoS with the phase transition is softer than the EoS of pure 
hadron gas, in some regions of thermodynamic parameters the EoS-PT has sound velocities 
of the same order or even larger than those in the EoS-HG. 

Below we consider collisions of gold nuclei. The hydrodynamic simulation is started at 
the stage when two cold Lorentz-contracted nuclei approach each other. The initial baryon 
and energy density profiles of each nucleus in its rest frame are described by the Woods- 
Saxon distribution. Due to the stabilizing effect of the Skyrme-like mean field, the initial 
nuclei may stay in equilibrium with vacuum at P = . They propagate without distortion 
until their density distributions essentially overlap. 

Below the beam axis is denoted by z and the x-axis is chosen along the impact parameter 
vector b. In our simulations we use linear interpolation of tables P{n, e) prepared with fixed 
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steps in e and n. Analogous tables for temperature (T), baryon (/i) and strange (/ig) 
chemical potentials are used for calculating hadronic spectra. The numerical solution of 
fluid-dynamical equations is obtained by using the relativistic version of the flux-corrected 
transport algorithm 

To calculate hadronic momentum distributions we apply the approximation of instan- 
taneous freeze-out: it is assumed that a sudden transition from the local equilibrium 
to collisionless propagation of particles takes place at some isochronous hypersurface 
{t = tfr = const). Within this approximation one can use the following formula jsl for the 
invariant momentum distribution of the hadronic species i 

where p'^ = {E, pY is the 4-momentum of the particle, y and Pt are, respectively, its longi- 
tudinal rapidity and transverse momentum, gi is the statistical weight of i-th hadrons. Plus 
and minus in Eq. (jl]) correspond, respectively, to fermions and bosons. Using conditions of 
chemical equilibrium one can express the particle's chemical potential Hi through the baryon 
and strange chemical potentials as follows 

/Xi = Bill + Sills , (5) 

where Bi and Si are, respectively, the baryon and strangeness number of species i . 

In addition to contributions of "thermal" nucleons and pions, which are calculated di- 
rectly by using Eqs. (4)-(5), we also take into account resonance decays, e.g. A — )• A^vr, 
p — 7- 27r,. . . We use the standard formulae for two-body decays of resonances in the zero width 
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approximation (the decays into more than two hadrons are treated approximately 
known resonances with masses up to 2 GeV are taken into account. 

The transverse collective flows of matter created in heavy-ion collisions, are rather sensi- 



tive to the EoS. Especially useful is the elliptic flow characterized by the parameter V2 |12 |. 
To discuss qualitatively possible differences between the EoS-PT and EoS-HG, below we 
calculate the so-called momentum anisotropy parameter ep which characterizes the flow 



asymmetry in the azimuthal plane. Following Ref. 



15| we define this quantity as 



I dxdy (T^^ - T^y) 
= -r , (6) 

dxdy (T^^ + Tyy) 
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where T^^^T^^ are the components of the energy-momentum tensor in the azimuthal plane 



0. The approximate relation ~ 2^2 has been obtained in Ref. 



13| from (2+1)- 



hydrodynamical simulations of Pb+Pb collisions at SPS and RHIC energies. 



III. RESULTS 



First, let us discuss thermodynamic characteristics of matter created in relativistic 
Au+Au collisions. Our goal is to find differences between the results obtained with EoS-HG 
and EoS-PT. Figures 1-2 represent dynamical trajectories of matter in the central box 
|2;|/7o < Ifm) produced in central (b=0) Au+Au collisions at different bombard- 
ing energies -Eiab- Shown are the results for n — e and /i — T planes. At -Eiab ^ 5 AGeV 
the calculation with EoS-PT predicts larger values of e and n as compared to EoS-HG. In 
this case one can also see longer life-times of states with maximal compression and delay 
in transition to expansion stages. Note that calculations with EoS-PT predict a zigzag-like 
behavior of the trajectories in the — T plane with slightly raising temperature in the MP 
as a function of time. It is important to note that the final states of hadronic matter in 
the central cell are very similar for two EoS used in the calculations. In other words, all 
differences in the early dynamics are practically washed out in the final stage. 

According to our analysis, especially interesting is the region of bombarding energies 
around E\ah ~ 10 AGeV. At such energies we find an enhanced sensitivity of collective flows 
and particle spectra to the phase transition. Figure [3] represents the time evolution of the 
momentum anisotropy Ep in semicentral Au+Au collisions at different E\^\^ . One can see that 
the momentum anisotropy is rather sensitive to EoS at i?iab ^ 20 AGeV. The asymptotic 
values of at -Eiab ^ 10 AGeV are larger in calculations with the deconfinement phase 
transition. This agrees with the similar conclusion of Ref. [l^ . 

To estimate sensitivity of these results to the choice of the freeze-out time, we determine 
the time moment when the energy density in a central box becomes smaller than a certain 
freeze-out value £fr. In Fig. [3] we mark the values of corresponding to different values of 
Sfj.. Figure 4 shows our results for the excitation function ep(i?iab)- The lines connect the 
values ep taken at freeze-out times corresponding to = 0.4 GeV/fm'^. In the case of EoS- 
PT we predict a non-monotonic dependence of ep(i?iab) with maximum at E'lab — 10 AGeV. 
On the other hand, the existing experimental data on the elliptic flow in Au+Au and Pb+Pb 
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collisions at AGS and low SPS energies do not show a non-monotonic behavior. This may 
imply that elliptic flows are strongly suppressed by viscosity effects in the considered energy 
domain. 

In Figs. 5-6 we present our results for proton and pion rapidity distributions in central 
Au+Au collisions at i^iab ~ 10 A GeV. The proton and vr^ distributions are obtained from 
nucleon and pion spectra by introducing the additional factors 1/2 and 1/3 respectively. 
The hadronic spectra were calculated by using Eqs. (4)-(5). The parameter tfj. is chosen to 
achieve the best flt of experimental data. To take into account the mean-field potential, we 
shift the chemical potentials of baryons by —U (n). In the case of EoS-PT we also introduce 
the excluded volume corrections. This is done by replacing /Xj — )■ /ij — VsPr, where Pk is 
the kinetic part of pressure sl]. According to our calculations, these corrections are more 
important for pion spectra. 

Figure |5] shows the rapidity distributions of protons in central Au+Au collision at 
-Eiab = 10.7 AGeV. The calculation with the phase transition predicts noticeably broader 
rapidity distributions. In this case the agreement with experimental data is better as com- 
pared with the EoS-HG. 

The IT" rapidity distribution for the same reaction is shown in Fig. [61 Again, calculations 
with the EoS-PT give broader distributions than those with the EoS-HG, however this 
difference is not so large as for protons. In the case of EoS-PT the observed pion spectra 
may be reproduced only if one takes a smaller freeze-out time as compared to protons. We 
expect that this differences would decrease if we would take smaller values of the excluded 
volume for pions. 



IV. CONCLUSIONS 



We have studied the sensitivity of the elliptic flow and particle spectra to the decon- 
flnement phase transition in heavy-ion collisions. Our analysis shows that at FAIR energies 
maximal values of energy- and baryon densities in the central box of the colliding system are 
signiflcantly larger if the QGP is formed at some intermediate stage of a heavy-ion collision. 
It it shown that the collective flow parameters are especially sensitive to the EoS around 
E'lab — 10 AGeV. At such energies the calculations with the deconflnement phase transition 
predict enhanced elliptic flows and broader proton rapidity distributions as compared with 
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the purely hadronic scenario. 
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FIG. 1: Time evolution of matter in central Au+Au collision at different bombarding energies -Eiab- 
Shown are values of energy and baryon densities averaged over the central box. Dashed and solid 
lines correspond to EoS-HG and EoS-PT, respectively. Numbers in circles and triangles give the 
cm. time in fm/c. Shading shows the MP region of deconfinement phase transition. 
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FIG. 2: Time evolution of matter (central cell) in the — T plane. Left and right panels correspond 
to central Au+Au collisions at £^iab = 10 and 20 AGeV, respectively. Dashed (solid) lines show the 
results for EoS-HG (EoS-PT). Dotted lines show the MP region. Numbers in circles and triangles 
give time values in fm/c. 
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FIG. 3: Time dependence of the momentum anisotropy in Au+Au collisions for different bombard- 
ing energies (6 = 7 fm). The solid (dashed) lines correspond to the EoS-PT (EoS-HG). Markers 
show time moments when energy density in the central box becomes lower than certain values efr 
indicated in the the right panel. 
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FIG. 4: Excitation function of momentum anisotropy in Au+Au collisions with b = 7 fm. The solid 
(dashed) lines correspond to the EoS~PT (EoS-HG). Arrows show possible shifts of Cp-values for 
different choices of freeze-out energy density between 0.2 and 0.6 GeV/fm^. 



1 1 1 1 1 1 


1 1 1 1 1 1 1 1 

Au+Au, b=7fm 












EoS-PT 




— EoS-HG 


1 , , , , 1 


efr=0.6 — efr=0.2 
, 1 , , , , 1 



11 




FIG. 5: Proton rapidity distributions for central Au+Au collision at -Eiab = 10-7 AGeV. Full sym- 



bols are experimental data 
midrapidity. 



16l4l8l]. Open symbols are obtained by reflection with respect to 




